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Abstract. We present an analysis of new Suzaku data and archival data from XMM-Newton of the cluster of galaxies Sersic 
159-03, which has a strong soft X-ray excess emission component. The Suzaku observation confirms the presence of the soft 
excess emission, but it does not confirm the presence of redshifted O vn lines in the cluster. Radial profiles and 2D maps derived 
from XMM-Newton observations show that the soft excess emission has a strong peak at the position of the central cD galaxy 
and the maps do not show any significant azimuthal variations. Although the soft excess emission can be fitted equally well 
with both thermal and non-thermal models, its spatial distribution is neither consistent with the models of intercluster warm-hot 
filaments, nor with models of clumpy warm intracluster gas associated with infalling groups. Using the data obtained by the 
XMM-Newton Reflection Grating Spectrometers we do not confirm the presence of the warm gas in the cluster centre with 
the expected properties assuming the soft excess is of thermal origin. The observed properties of the soft excess emission 
are consistent with the non-thermal interpretation. While the high density of relativistic electrons associated with the peak of 
the soft emission in the cluster centre might have been provided by an active galactic nucleus in the central cD galaxy, the 
underlying population might have been accelerated in diffuse shocks. 
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1. Introduction 

In the past 10 years, extreme ultraviolet (EUV) and soft X-ray 
observations (with EUVE, ROSAT, BeppoSAX, XMM-Newton) 
unveiled emission in excess of that expected from the ther- 
mal intracluster medium (ICM) in a number of clusters of 
galaxies. The first papers reporting the discovery of the ex- 
cess EUV emission from clusters interpre ted their finding as 



thermal emission from warm diffuse gas dLieu et al. I ll996albt 



Mitt az et al. 19981) . However, the large mass of the warm gas 
required to explain the soft excess emission, which would 
cool very rapidly at the derived tem peratures, was a prob- 



lem for the model ( Mittaz et al. 1998b . As an alternative ex 



planation for the cluster soft excess inverse Compton (IC) 
emission by cosmic-ray electrons scattering off the cosmic 
microwave background (CMB) was proposed dHwand["l997l: 



Sarazin & Lieu 



19987 



Lieu et al 



Ensslin & Biermann] [l998l 

1999). Kaastra et al.l d2003l) reported the discovery of soft ex- 
cess and Ovn line emission around five clusters observed 
with XMM-Newton. They attributed this component to emis- 
sion from intercluster filaments of the Warm-Hot Intergalactic 
Medium in the vicinity of the clusters. [Nevalainen et alj Q2006) 



reanalized the XMM-Newton data of 4 clusters with a reported 
soft excess emission using the newest calibration and found 
that the XMM-Newton EPIC instruments differ on the mag- 
nitude of the soft excess in all clusters. While EPIC/MOS 
data still show soft exc ess emission in all clusters reported 
by iKaastra et al.l ([2003 ), EPIC/pn data show the soft excess 
emission only in the cluster Sersic 159-03 and the redshifted 
Ovn line emission is n ot observed in the reprocessed data. 
Nevalainen et al. I d2006h conclude that the possibility that the 
reported O vn line emission in the outskirts of these clusters is 
due to heliosphe ric or geocoronal charge exchang e emission, 
as suggested by iBregman & Llovd-Daviesl d2006l) . cannot be 
ruled out. 

The rich southern cluster of galaxies Sersic 159-0 3, also 
known as ACO S 1 1 01, was discover ed by [Sersic dl974l) . It is a 
nearby (z=0.0564: lMaia et al.lfl987h X-ray bright cluster with 
a lum inosity of L x = 5.35 X 10 44 erg s 1 in the 0.5-2.0 keV 
band ddeGrandi etalJll999h . The cluster is relaxed with no 



obvious peculiarities in the X-ray morphology dKaastra et al 
200 II) . The temperature profile peaks at kT = 2.7 keV at a 



Send offprint requests to: N. Werner, email n.werner@sron.nl 



radius of ~ 2.5' from the core dde Plaa et al. I l2006l) . The tem- 
perature drop in the core is relatively modest and the temper- 
ature outside the radius of 4' drops rapidly by at least a factor 
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of ~ 2 dKaastra et al.ll200Ulde Plaa et al.ll2006l) . Sersic 159-03 
has the strongest soft excess emission of all clusters observed 
with XMM-Newton, with the detection well above the calibra- 
tion uncertain ty. Its soft excess was di scovered using ROSAT 
PSPC data by [Bonamente et alj (1200 lb an d confirmed by two 



independent XM M-Newton observations (IKaastra et al.l 12003 
Bonamente et al. 2005 ; Nevalainen et al. 20061). Analysing a 



deep (122 ks) XMM-Newton observation, Bon amente et al 



(2005) show that strong soft excess emission is detected out 
to the radial distance of 12' (0.9 Mpc). They interpret their 
data using two models: by invoking a warm reservoir of ther- 
mal gas, and by relativistic electrons that are part of a cos- 
mic ray population. In their conclusions they slightly favor the 
ther mal interpretation. B ased on the analysis of the same data 



set, de Plaa et alj |2006) propose as a potential source of the 



soft excess non-thermal emission arising from IC scattering be- 
tween CMB photons and relativistic electrons, that are acceler- 
ated in bow shocks associated with ram pressure stripping of 
in-falling galaxies. 

Here we analy se the soft excess in Sersic 159-03 using a 
new deep Suzaku dMitsuda et al. 20071) observation , together 
with two archival XMM-Newton (Jansen et al. 2001) observa- 
tions. We analyse the Suzaku XIS1 (IKoyama et al.ll2007l) data 
in order to obtain a confirmation of the soft excess with an inde- 
pendent instrument. We compare the observed soft excess flux 
in the three analysed observations and determine the systematic 
uncertainty on its value. Furthermore, Suzaku has a superior 
spectral redistribution function at low energies, which allows 
us to resolve the oxygen line emission with unprecedented ac- 
curacy. The detection of O vn line emission with Suzaku would 
prove the thermal origin of the soft excess emission and deter- 
mining its redshift would allow us to confirm its cluster origin. 
We also take advantage of the excellent statistics of the archival 
XMM-Newton observations, which allow us to analyse the spa- 
tial distribution of the soft excess emission in the cluster. 

The paper is organised in the following way. In section 2, 
we describe the Suzaku and XMM-Newton data of Sersic 159- 
03 and discuss the spectral modeling. In section 3, we fit the 
soft excess emission detected by Suzaku with thermal and non- 
thermal models, we compare the soft excess flux detected by 
Suzaku with the flux detected during two XMM-Newton ob- 
servations. Furthermore, we exploit the higher signal-to-noise 
and good spatial resolution of the XMM-Newton observations 
to investigate the spatial structure of the soft excess emission by 
analysing radial profiles and 2D maps, and we search for possi- 
ble lines from warm gas in the cluster core with the Reflection 
Grating Spectrometers (RGS) on XMM-Newton. In section 4, 
we discuss the physical interpretation of the results and fi- 
nally, in section 5, we summarise the main conclusions. In the 
Appendix we describe the details of the calibration of the con- 
taminating layer on the optical blocking filter of the XIS 1 de- 
tector on Suzaku. 

Throughout the paper we use Hq = 70 km s _1 Mpc -1 , Q.m = 
0.3, Qa = 0.7, which imply a linear scale of 73 kpc a rcmin 1 
at the cluster redshift of z = 0.0564 dMaia et al.lll987b . Unless 
specified otherwise, all errors are at the 68% confidence level 
for one interesting parameter (Ay 2 = 1). Upper limits are at the 




Fig.l. Suzaku XIS image of Sersic 159-03 in the 1-4 keV 
band. The data obtained by the four XIS detectors were com- 
bined. The white circles indicate the extraction regions for 
the spectral analyses: circular region with a radius of 3' (the 
smaller circle), and annulus with inner and outer radii of 3' and 
8', respectively. 



2cr confidence level. The elementa l abundances ar e given with 
respect to the proto-solar values of lLoddersI (120031) . 



2. Observations and data analysis 

2.1. Suzaku data 

Sersic 159-03 was observed with Suzaku between April 26-28, 
2006 with a total exposure time of 64 ks. 

We analyse the data obtained by the back illuminated XIS 1 
detector, using event files of version 1.2 products. We only use 
the XIS 1 data and not the remaining 3 front illuminated XIS 
detectors (0, 2, 3) in order to reduce the systematic uncertain- 
ties. The XIS1 detector has the largest effective area and is best 
calibrated at the low energies of all XIS instruments. XIS1 has 
also the best studied radial profile of the contamination of the 
optical blocking filter. 

The XIS CCD camera covers the 0.2-12 keV energy range 
with an energy resolution of 40 eV at 0.5 keV and 130 eV at 
5.9 keV. It has a square field of view of 18'xl8'. Its advantage is 
the absence of a low energy tail in the pulse-height distribution 
function, which makes it an excellent instrument to study line 
emission at low energies. 

The image of Sersic 159-03 in the keV band is shown 
in Fig. [T] Vignetting effects were not taken into account for this 
image and the background events were not subtracted. Suzaku 
XIS 1 has a relatively b road point spread funct ion with a half- 
power diameter of 2.3' (ISerlemitsos et al .120071) . Therefore, the 
extraction regions also have to be correspondingly large. The 
circles indicate two regions, that we used for spectral analysis. 
The spectra were extracted from a circular region with a radius 
of 3' centred on the cluster and from an annulus with inner and 
outer radii of 3' and 8', respectively. We extracted the spectra 
using XSELECT distributed with HEASOFT 6.1.1. 
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To subtract the instrumental background, we extracted 
spectra from night Earth observations. The level of the instru- 
mental background is anti-correlated with the cut-off rigidity 
(COR) value at the position of the satellite. We extracted night 
Earth spectra for each 1 GeV c _1 interval of COR and added 
them weighted by exposure time for the given COR range in 
the source observation. The background spectra were extracted 
from the same detector region as the source spectra in order to 
avoid possible systematic effects due to spatial variation of the 
instrumental background. Cosmic background components are 
included in the spectral fitting phase (see Sect. 12.41 ). 

The spectral redistribution files (RMF) were created us- 
ing 'xisrmfgen'. The ancillary response files (ARF) were 
created using 'xissi marfgen', which is based on ray-tracing 
(Ishi saki et al. 20071) . We created a separate ARF for both ex- 
traction regions. For the incident flux distribution, we adopted 
the XIS image of the cluster in the 1 .0^4.0 keV band. The ef- 
fective area of the XIS detectors below 1 keV is affected by 
carbon and oxygen contamination of the optical blocking filter. 
The contamination is time and position dependent. To model 
the column and the position dependence of the contaminant, 
we use the "ae_xil_contamiJ20061016.fits" contamination ta- 
ble, which contains the radial profiles and the best estimated 
values of the contaminating column for different time periods. 
For the precise value of the contaminating column density in 
the time of our observation, we use the values determined using 
simultaneous calibration observations of the blazar PKS 2155- 
304 performed by Chandra LETGS/HRC, XMM-Newton and 
Suzaku on May 1st 2006, only 3 days after our observation (see 
the Appendix for the details of the analysis). We found that 
the column densities of the contamination that are suggested 
by the XIS team (Nc = 3.5 x 10 18 ctrr 2 , N = 5.9 x 10 17 

cm 2 ) are not consistent with the calibration observation of 2.3. Spectral analysis 



while the former value is too low and is not consistent with the 
values obtained from the analysis of the cross-calibration data, 
the second value is only slightly higher than that obtained by 
the cross-calibration of EPIC/pn and XIS 1 . 

We inspected the solar proton flux measured with the ACE 
and Wind satellites during our observation. The proton flux in- 
creases in the last 24 ks of the observation to ~10 9 s _1 cm 



level a t which the XIS1 spectrum is affected in lFuiimoto et al 



2 , the 



d2007h . In order to investigate the possibility of the contami- 
nation of our data by charge exchange emission in the Earth's 
magnetosphere, we extracted a spectrum using only the first 40 
ks, with a low proton flux. We compared it with the spectrum 
extracted from the remaining 24 ks of the Suzaku observation. 
We verified that the spectrum did not change and that the ele- 
vated soft proton level in the second part of the observation did 
not cause a detectable contamination by charge exchange emis- 
sion. Hence we used the whole observation in our analyses. 

2.2. XMM-Newton data 

Sersic 159-03 was observed with XMM-Newton on May 5, 
2000 and on November 20-2 1 , 2002 with a total exposure time 
of 60 ks and 122 ks, respectively. The data were reduced with 
the 7.0.0 version of the XMM Science Analysis System (SAS). 
The data processing and the subtraction of the instrumental 
background for the European Pho ton Imaging Camer as (EPIC) 
was performed as described in Ide Plaa et alJ d2006h . Cosmic 
background components in the EPIC data are included in the 
spectral fitting phase. The R GS data were processed as de- 
scribed in ldePlaaetai1 ( l2006l) . 



PKS 2155-304. The carbon and oxygen contaminating columns 
determined using the simultaneous observation of Suzaku and 



Chandra LETGS/HRC are Nc = (3.78 ± 0.05) x 10 ls cm z and 
No = 6.30 + 0.08 x 10 17 ctrT 2 , for the centre of the field of 
view. Since the most probable composition of the contaminant 
is C24H38O4, the C/O ratio was fixed in the fitting process to 
6 (see Appendix for the details). The contaminating columns 
determined using the simultaneous Suzaku and XMM-Newton 
EPIC/pn observation are Nc = (4.11 ± 0.06) x 10 18 cirT 2 and 
No = (6.2 + 0.5) x 10 17 ctrT 2 , for the centre of the field of view. 
This analysis shows, that there is ~10% uncertainty in the value 
of the contaminating column. 

In order to provide a more conservative estimate for the 
soft excess, we will adopt in this paper the values for the con- 
taminating column in the centre of the field of view deter- 
mined by the simultaneous observation of Suzaku and Chandra 
LETGS/HRC: N c = 3.78 x 10 18 ctrT 2 and N = 6.30 x 
10 17 ctrr 2 . These values correspond to the nominal values in 
the standard contamination file for July 1st 2006. To illus- 
trate how sensitive the best fit soft excess flux is to the pre- 
cise value of the contamination, we will also compare our re- 
sults in Sect. 13.21 with those obtained assuming the nominal 
column densities in the contamination file for the date when 
our observation was performed [Nc = 3.5 x 10 18 ctrT 2 ) and 
for September 1st 2006 (7V C = 4.2 x 10 18 ctrT 2 ). We note that 



For the spectral analysis we use the SPEX package 
dKaastra et al. 19961) . We fix the Galactic absorption in our 
model to the value deduced from Hi data Ng = 1.79 x 
10 20 cm -2 dDickev & Lockmanll990l) . We use the lVerner eTal 



dl996 ) cross-sections in our model of Galactic absorption. For 
the spect ral fitting of the c luster spectra, we use the wdem 



model dKaastra et al.l |2004), which prov ed to be the most 
successful in fitting cluster cores (e.g . iKaastra et al 



dePlaaetal. 2005; Werner et al. 2006; dePlaaetal. 2006). 



Ka; 
Idel 



2004 



This model is a differential emission measure (DEM) model 
with a cut-off power-law distribution of emission measures ver- 
sus temperature. In our analysis, we fit both a thermal plasma 
model (MEKAL) and a non-thermal power-law model to the 
soft excess emission. In our spectral fits both the cluster and 
the soft excess emission are free parameters. 

Unless specified otherwise, we use the 0.4-7 keV energy 
range for the spectral analysis of the Suzaku XIS 1 data and the 
0.4-10 keV band for the XMM-Newton data. The Suzaku and 
the two XMM-Newton data sets are fitted separately. 

2.4. Modeling of the X-ray background emission 

We correct for the Cosmic X-ray Backgro und (CXB) during 
spectral fitting. iKuntz & Snowdenl ((2000) distinguish 4 dif- 
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Table 1. The unabsorbed soft Galactic foreground flux in units 
of 1CT 12 erg s _1 cm~ 2 deg" 2 in the 0.3-10 keV band. The col- 
umn RA SS shows the fluxes deriv ed from the parameters re- 
ported by Kuntz & Snowden ( 2000l) . In columns Suzaku A and 
B, we rep ort the value s deriv ed from two offset pointings near 
A2218 bv lTakeietaHd2007l) . The values in column S159-03 
were derived using XMM-Newton data in the 9'-12' a nnulus 
around the cluster Sersic 159-03 bv lde Plaa et al.1 d2006l) . 



Component 


RASS 


Suzaku A 


Suzaku B 


S 159-03 


LHB kr 


0.082 


0.08 


0.08 


0.07 


LHB flux 


1.22 


1.74 


3.35 


2.23 


SDC kr 


0.068 








SDC flux 


0.87 








HDC kr 


0.127 


0.16 


0.25 


0.20 


HDC flux 


2.44 


1.85 


2.06 


12.20 


total soft flux 


4.53 


3.59 


5.41 


14.43 



ferent background/foreground components: the extragalactic 
power-law (EPL), the local hot bubble (LHB), the soft dis- 
tant component (SDC) and the hard distant component (HDC). 
The EPL component is the integrated emission of faint dis- 
crete sources, mainly distant Active Galactic Nuclei (AGNs). 
The LHB is a local supernova remnant, in which our Solar 
System resides. It produces virtually unabsorbed emission at 
a temperature of ~10 6 K. The soft and hard distant components 
originate at larger distances. They might be identified with the 
Galactic halo, Galactic corona or the Local group emission and 
are absorbed by almost the full Galactic column density. In the 
first column of Table Q] we show the fluxes and temperatures 
for the G alactic foreground derived from the parameters re- 
ported by iKuntz & Snowdenl (2000), based on the Rosat All 
Sky Survey data. Using the spectral band above 0.4 keV we 
can not reliably distinguish the SDC emission from the LHB 
component. Therefore, at temperatures below 0. 1 keV we only 
consider the contribution of the LHB. 

Using data obt ained by XMM-Newton, 
De Luca & Molendil d2004l) found that the photon index 
of the EPL is T = 1.41 + 0.06 and its 2-10 keV flux is 



(2.24 + 0.16) x lO" 11 erg ctrT 2 s~' deg~ 2 (90% confidence 
level). By fitting the spectra extracted from the 9'— 12' re- 
gion arou nd Sersic 159-03 and assuming a photon index of 
T = 1.41. ldePlaaet~ai1d2006l) found a 2-10 keV EPL flux of 
2.26 x 10~ n erg ctrT 2 s deg~ 2 (corresponding to a 0.3-10.0 
keV flux of 3.14 x 10~ n ergs cm~ 2 s _1 deg -2 ). We adopt 
these values in our model of the EPL background. While the 
measured mean variation in the EPL emission in the 2-10 keV 
band in the sky is only ~ 3.5%, the intensity of the soft 
emissio n in the 0.2-1 .0 keV band varies by ~ 35% from field 
to field dLumb et al.ll2002i: iKaastra et alj|2003l) . Unfortunately, 
the emission from Sersic 159-03 fills the entire field of view 
of both Suzaku and XMM-Newton, which does not allow us 
to determine the properties of the local background directly 
from an annulus in the outer part of the field of view of the 
instruments. As shown in Table PT1 Ide Plaa etaL ( 20061) found, 
by fitting the spectra extracted from the 9'— 12' region around 
the cluster, a flux of the HDC emission which was more than 



5 times highe r than the HDC flux in the soft band measured 
by ROSAT dKuntz & Snowdenl I200C ), and determine d from 
two Suzaku blank field observations ( Takei et al. I l2007h . Since 
the flux at low energies in the outer parts of the cluster is not 
consistent with the values determined in other fields, large 
part of it may still be associated with the soft excess emission, 
which we are looking for. 

We inspected the ROSAT all sky survey data in the 3/4 keV 
band (0.5-0.9 keV; which is less affected by the variations in 
the Galactic absorption) and in the 1/4 keV band (0. 1-0.4 keV) 
around Sersic 159-03 and found that the variations in the soft 
emission around the cluster are small and there is no indication 
of a strongly elevated level of Galactic foreground. The ROSAT 
count rates in the vicinity of Sersic 159-03 are consistent with 
those in the Suzak u offset pointing B at the cluster Abell 2218 
dTakei et al.l 120071) . Therefore, for the soft foreground emis- 
sion, we adop t the fl uxes reported for the offset pointing B 
by iTakei et al. (2007). The sum of the adopted fluxes of the 
soft foreground components is ~ 20% high er than the average 
soft X -ray background flux determined by IKuntz & Snow den 
(2000) from the ROSAT observations. 



3. Results 

3.1. Thermal and non-thermal models of soft excess 
emission 

We analyse the Suzaku XIS 1 spectra extracted from the two re- 
gions shown in Fig. 1 : a circular region with a radius of 3', and 
annulus with inner and outer radii of 3' and 8', respectively. We 
fit the cluster emission with the multitemperature wdem model 
(see Sect. 12.3b . The excess emission at low energies can be for- 
mally fitted as strongly subgalactic absorption column density 
(A^h). When fitting the cluster spectra with free A^h, we obtain 
for the inner extraction region an absorption column density 
of iV H = (1.15 + 0.20) x 10 2() ctrT 2 and for the outer region 
we obtain a 2<x upper limit of A^h < 1.4 x 10 19 crrT 2 . Both 
values are well below the Galactic value toward the cluster 
A^h = L79 x 10 20 cm" 2 (for an extensive discussion on ab- 
sorption toward Sersic 159-03 see iBonamente et aTlbOOlb . In 
the rest of the analysis, we fix the A^h in our spectral model to 
the Galactic value. 

In order to describe the soft excess emission, we try to fit 
both thermal and non-thermal models. Both the parameters of 
the cluster component (ICM emission measure, temperature 
structure, and abundances) and of the soft excess emission are 
free to vary in the fitting process. When fitting the soft excess 
emission with the thermal model, we fix the abundances of the 
warm-hot gas other than oxygen to 0.3 solar and fit its emis- 
sion measure, temperature and oxygen abundance. In the non- 
thermal model, the power-law emission caused by relativistic 
electrons is assumed to contribute both to the soft excess and 
some part of the higher energy emission. When fitting the soft 
excess emission with a non-thermal component, we fit the pho- 
ton index and the normalisation of the power-law. 

In Table |2] and in Fig. [2] we show the best fit results using 
both thermal and non-thermal models to fit the soft excess. We 
also show the best fit results for the cluster emission, including 
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0.5 1 2 5 0.5 1 2 5 

Energy (keV) Energy (keV) 

Fig. 2. Left panel: Suzaku XIS 1 spectrum of the 3'-8' annulus around the cluster. The total model and its components are plotted 
separately. The soft excess emission is fitted with a thermal model. The dash-dotted line indicates the subtracted instrumental 
background, the dotted line indicates the cosmic X-ray background, the dashed line shows the thermal emission associated with 
the soft excess, the dash-dotted-dotted-dotted line indicates the cluster emission and the full line indicates the total model. Right 
panel: The same spectrum as on the left panel, with the soft excess modelled as non-thermal emission. 



the best fit abundance values for 6 elements. When fitting the 
data from the inner extraction region the oxygen abundance of 
the soft excess component is not constrained by the data and we 
fix its value to 0.3 solar together with the other abundances. The 
best fit temperature of the warm component does not depend on 
the adopted abundance values. For lower adopted metallicity, 
the best fit emission measure of the warm gas will increase and 
for higher metallicities it will decrease. The soft excess flux in 
the cluster centre is higher than in the outer extraction region, 
which indicates that the soft excess peaks at the cluster core. 

As we show in Table [2] and in Fig. [2] thermal and non- 
thermal models fit the data statistically equally well. We note 
that the reduced \ 2 values of our best fit models are rela- 
tively high because we did not include systematic errors on 
the data in the spectral fitting process and at the high number 
of counts of our dataset the systematic errors dominate over 
the statistical errors. Clear detection of a redshifted O vn line 
emission would be an unambiguous proof for the thermal ori- 
gin of the soft excess emission. Since XIS1 has a much bet- 
ter spectral redistribution function at low energies than EPIC, 
we use its data for the outer extraction region to search for 
the redshifted O vn line emission. In order to model the broad 
soft excess we fit a power-law, and subsequently add a nar- 
row Gaussian at the expected energy of the O vn intercombi- 
nation line at the cluster redshift. We find that the power-law 
describes the soft excess emission sufficiently well and there is 
no need for the additional Gaussian. This way, we find a for- 
mal upper limit on the O vn line emission from the cluster of 
1.69 x 1CT 7 pho tons cm" 2 s" 1 arc min 2 , which is comparable 



to that found by lTakei et all (120071) for the cluster Abell 2218. 
This demonstrates the sensitivity of XIS for lines at low ener- 
gies. Although the thermal model of the excess emission which 
already contains the Ovu line emission (with a line flux of 



XI 

o 



XI 

o 



XI 
> 



O 



Suzaku XIS1 
O EPIC/MOS 
EPIC/pn 




in 
d 
I 



0.5 



Energy (keV) 



(6.4 + l)xl0 photons cm s arcmin ) is consistent with 
the data, we can not confirm the presence of the redshifted O vn 



Fig. 3. The residuals of the best fit cluster model in the 0.9-7.0 
keV band to the spectra obtained by Suzaku XIS 1 , EPIC/MOS, 
and EPIC/pn from the 3'-8' annulus around the cluster, extrap- 
olated to low energies. All instruments show a soft excess emis- 
sion below 1 keV. 

line emission which would be the clear proof for the thermal 
origin of the soft excess. In subsection 13.51 we show that the 
O vii line emission is also not present in the RGS data, as would 
be expected given the spatial distribution of the excess emission 
(see Sect. EDI 

3.2. Comparison with XMM-Newton soft excess 
detections and the systematic uncertainties 

In order to determine the level of consistency between XMM- 
Newton EPIC and Suzaku XIS 1 below 1 ke V, we fit separately 
for each instrument the cluster spectra extracted from the 3'- 
8' annulus in the 0.9-7.0 keV band with the multitemperature 
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Table 2. Fit results for the Suzaku XIS1 data of Sersic 159-03. 
We fit the multi-temperature wdem model to the cluster emis- 
sion, combined with both thermal and non-thermal models to 
describe the soft excess in a circular region with a radius of 3' 
centred on the cluster core and in an annulus with inner and 
outer radii of 3' and 8', respectively. For the Galactic absorp- 
tion we use a value of 1.79 x 10 20 crrT 2 . Emission measures 
(Y = f n e riudV) are given in 10 66 crrT 3 and the power-law 
normalisations (A^nt) are given in 10 51 photons s" 1 keV" 1 . The 
absorbed fluxes (F) were determined in the 0.3-10.0 keV band 
using an A^h value of 1.79 x 10 20 ctrT 2 , and they are given in 



10 ergs s cm deg Abundances ar e given with respect 
to the proto-solar values of lLoddersI (12003I) . 



Par. 



F a 

kT 

IY1 me; 

o 

Si 

s 

Ar 
Ca 
Fc 
Y 

1 exc 

kT 
O 

JVnt 

r 



0'-3' 
therm. 



0'-3' 
non-therm. 



3'-8' 
therm. 



3'-8' 
non-therm. 



15.8 ±0.2 
30.68 ± 0.38 
2.61 ±0.08 
0.46 ±0.12 
0.48 ± 0.04 
0.38 ± 0.05 
0.31 ±0.13 
0.72 ±0.18 
0.58 ± 0.02 

0.6 ± 0.2 
0.31 ±0.05 
0.3 



14.9 ± 0.5 
29.03 ± 1.02 
2.56 ±0.09 
0.64 ± 0.08 
0.50 ± 0.04 
0.41 ± 0.06 
0.35 ±0.14 
0.78 ±0.19 
0.60 ± 0.02 



3± 1 

2 11 +0 ' 27 



3.71 ±0.09 
1.03 ±0.03 
2.41 ±0.14 

< 0.30 
0.21 ±0.07 

< 0.52 

< 0.32 

< 0.62 
0.39 ± 0.03 
1.25 ±0.25 
0.19 ±0.01 
0.11 ±0.04 



3.2 ±0.3 
0.92 ± 0.06 
2.37 ±0.19 
0.51 ±0.15 
0.25 ± 0.08 
0.15 ±0.11 

< 0.36 

< 0.46 
0.41 ± 0.04 



2.0 ±0.7 
2 41+ - 33 



0.48 ±0.13 2.42 ±0.76 0.08 ±0.01 0.23 ± 0.08 



203/128 



212/128 



154/108 



164/108 



Table 3. The absorbed soft excess flux in units of 10~ n erg s" 1 
cm -2 deg -2 in the 0.3-1.0 keV band, determined using an /Vh 
value of 1.79 x 10 20 crrT 2 . The soft excess was modelled with 
a thermal model. We show the soft excess in two extraction 
regions. The datasets obtained during the first and second ob- 
servation with XMM-Newton are labelled EPIC 1 and EPIC 2, 
respectively. For the XIS 1 instrument, we report the soft excess 
fluxes obtained for 3 different contaminating columns. See the 
text for the details. 



Dataset 


0' 


-3' 


3'- 


8' 


EPIC 1 


5.95 


±0.35 


0.86 ± 


0.08 


EPIC 2 


6.03 


±0.23 


1.12 ± 


0.06 


XIS A 







0.23 ± 


0.07 


XIS B 


4.8 


± 1.3 


0.75 ± 


0.08 


XIS C 


7.4 


± 1.3 


1.29 ± 


0.10 



wdem model (see Sect. I2.3l l and plot the residuals for each in- 
strument extrapolated to low energies. The residuals for Suzaku 
XIS 1 , EPIC/MOS, and EPIC/pn are shown in Fig.|3] We clearly 
see that all instruments consistently show a soft excess emis- 
sion below 1 keV. We note that the observed excess is well 
above the calibration uncertainties. Observations of calibration 



sources with XMM-Newton do not show significant positive 
residuals below 1 keV. 



In order to determine the flux in excess of the tail of the 
hot bremsstrahlung emission below 1 keV, we fit the soft ex- 
cess emission with a thermal model (kT = 0.19 keV, metal- 
licity of 0.3 solar and emission measure as a free parameter), 
as a convenient empirical model. In order to provide a consis- 
tent comparison of the soft excess level between the different 
observations and between the two extraction regions we fit the 
same model to all datasets. The ICM emission measure, tem- 
perature structure and abunda nces are free parameter s in the 
fitting process. As repor ted by lBonamente et al.l (120051) and by 
Nevalainen et al. (2006), there is a disagreement between the 
XMM-Newton instruments on the level of the soft excess. We 
find that the spread in the best fit soft excess flux determined 
by fitting separately the data obtained by the individual instru- 
ments during the two XMM-Newton observations is smaller 
than 30%. In Table [3] we report the intensity of the soft ex- 
cess in the 0.3-1.0 keV band, obtained from individual obser- 
vations by XMM-Newton and Suzaku, in the two extraction 
regions shown in Fig. 1. The two observations with XMM- 
Newton EPIC are indicated as "EPIC 1" and "EPIC 2", respec- 
tively. In order to indicate the possible systematic uncertainties 
in the results obtained with XIS1, for this instrument we report 
the soft excess fluxes obtained for 3 different contaminating 
columns: for the nominal value in the time of our observation 
(indicated as "XIS A"; N c = 3.5 x 10 18 ctrT 2 ), for the best fit 
contaminating column indicated by the simultaneous calibra- 
tion observation with Suzaku and Chandra LETGS/HRC (indi- 
cated as "XIS B"; N c = 3.78 x 10 18 cm" 2 ), and for the upper 
limit of the contaminating column indicated by the calibration 
observations (indicated as "XIS C"; N c = 4.2 x 10 18 cm -2 ). 
We see that all observations show the soft excess emission in 
the 3'-8' annulus, and except "XIS A", the 0.3-1.0 keV fluxes 
are » 1.0 x 10~ n erg s" 1 deg -2 . Except of "XIS A", the soft 
excess flux in the cluster centre is a factor of ~6 higher than 
in the outer extraction region. For "XIS A", we detect in the 
cluster core excess absorption (A^h = 2.25 x 10 20 cirT 2 ) instead 



of excess emission. However, as we described in the Sect. 12.11 
and as we show in the Appendix, the simultaneous calibration 
observation with Suzaku, Chandra LETGS/HRC and XMM- 
Newton performed only 3 days after our observation allowed 
us to rule out the low value for the contaminating column with 
which the "XIS A" values were determined. The analysis of 
the simultaneous calibration observations shows that the con- 
tamination is between the values that were assumed in the de- 
termination of the fluxes indicated as "XIS B" and "XIS C". 
Taking into account the differences in the soft excess flux below 
1 keV shown in Table [3] we conclude that the soft excess flux 

- 2 deg -2 and(1.0±0.3)xl0- n 
deg 2 in the inner and outer extraction region, re- 
spectively. The ~20% uncertainty in the Galactic foreground 
flux introduces an additional ~10% systematic uncertainty to 
the soft excess flux determined for the 3'-8' extraction region. 



is (6.0 + 2.0) x 10" 11 erg s" 1 cm 
erg s" 1 cm" 2 
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Fig. 4. Left panel: Radial profiles of the flux of the hot ICM and of the soft excess emission in the 0.3-1.0 keV band derived 
from XMM-Newton data. The soft excess is modelled as thermal emission (kT = 0.19 keV). Central panel: Radial profile of 
the fraction of the soft excess flux relative to the hot ICM flux in the 0.3-1.0 keV band. The soft excess is modelled as thermal 
emission. Right panel: The radial distribution of the deprojected electron density of the warm gas, assuming the soft excess is of 
thermal origin and assuming spherical symmetry of its distribution. 



3.3. Radial distribution of the soft emission 

We determine the radial profile of the soft component using 
XMM-Newton EPIC spectra extracted from circular annuli, 
centred on the core of the cluster. We use annuli with outer 
radii of 0.5', 1', 2', 3', 4' and 6'. The cluster emission in each 
annulus is modelled with the multi-temperature wdem model 
(see Sect. EH). 

To account for projection effects, we deproject our spectra 
under the assumption of spherical symmetry. First, we fit the 
spectrum in the outermost (6th) annulus and determine the nor- 
malisation, temperature and iron abundance of the hot cluster 
gas and the normalisation of the soft excess component. We de- 
termine which fraction of this emission is projected in front and 
behind of the 5th annulus. We include the projected emission 
of the 6th shell, with fixed parameters, in the model for the 5th 
annulus. In the same way, in the model for each annulus, we in- 
clude the projected hot ICM and soft excess emission from all 
intervening shells. For the soft excess component, we try to fit 
both a power-law and a thermal plasma model. For the temper- 
ature of the thermal model and for the differential photon index 
of the non-thermal model, we adopt the best fit values obtained 
by fitting the spectrum extracted from the 3 '-8' region around 
the cluster (kT = 0.19 keV and F = 2.41), where the relative 
contribution of the soft excess emission is high and allows for 
the best determination of its spectral properties. 

On the left panel of Fig. [4] we show the radial profile of the 
projected soft excess and hot ICM flux in the 0.3-1 .0 keV band. 
We model the soft excess emission with a thermal model. On 
the central panel we show the radial profile of the fraction of 
the soft excess flux relative to the ICM flux in the 0.3-1.0 keV 
band. Fig.|4]shows that while the soft excess emission peaks at 
the core of the cluster, the soft excess flux relative to the cluster 
emission is increasing with radius beyond ~3'. 

Using the thermal description of the soft excess emission 
and assuming that the warm gas resides in the same shells as 
the hot ICM plasma, we calculate the electron density of the 
warm gas in each shell (see the right panel in Fig.[4]l. Using this 
deprojected thermal electron density profile, we obtain a total 
mass for the warm gas in the central region within a radius of 



2' (150 kpc) of M warm = 4.4 x 10 11 M . The mass of the hot 
ICM in the same region is M hot = 2.6 x 10 12 M G . 

On the left panel of Fig.|5]we show the relative contribution 
of the non-thermal flux between 0.3-10 keV, assuming the soft 
excess is of non-thermal origin. In order to illustrate the uncer- 
tainties in the radial profiles, we indicate the values obtained by 
EPIC/MOS and EPIC/pn separately. The plot shows that there 
are systematic differences in the best fit flux of the soft excess 
component determined using EPIC/MOS and EPIC/pn, how- 
ever, both instruments show the same radial trends. The best fit 
values of the soft excess in the radial bins have a systematic 
uncertainty of the order of ~40%. The secure determination of 
the soft excess flux depends critically on the instrument cali- 
bration. 

Using the non-thermal description for the soft excess emis- 
sion, we calculate the average density of the hot ICM (see the 
central panel in Fig. [5]) and the luminosity of the power-law 
component in each of the investigated shells. We use the depro- 
jected luminosities of the power-law component to determine 
the energy in the populatio n of relativistic electrons for each of 
the shells JLieu et alJll999h : 

£ f = 8x 10 61 L 42 1 y r" y r ergs, (1) 

where L42 is the luminosity of the non-thermal radiation in 
units of 10 42 ergs s , y m in and y max are the Lorentz factors cor- 
responding to the lower and higher limits of the energy range 
of XMM-Newton (0.3-10 keV) and p is the index of the differ- 
ential number electron distribution. For the Lorentz factors we 
adopt values y m j n = 600 and y max = 3500, using the relation 
y = 300(£/75 eV) 1/2 , where E is the energy to which a CMB 
photon is up-scattered by a relativistic electron with a Lorentz 
factor of y. The relationship between the differential photon in- 
dex F of the observed power-law emission and the index of the 
electron distribution isp — -1 +2F. We calculate the density of 
the relativistic electrons in the given shell by dividing the en- 
ergy density in relativistic electrons by <y>m e c 2 , where <y> 
is the average Lorentz factor (« 1000) and m e is the rest mass 
of an electron. 
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Fig. 5. Left panet.Radial profile of the relative fraction of non-thermal emission in the 0.3-10 keV band, if the soft excess is 
modelled as non-thermal emission. In order to illustrate the systematic uncertainties, we plot the EPIC/MOS and EPIC/pn data 
separately. Central panel: The radial distribution of the deprojected electron density of the hot ICM plasma. The soft excess 
was fitted with a power-law model. Right panel: The radial distribution of the deprojected density of the relativistic non-thermal 
electrons. 
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Fig. 6. Left panel: Radial profile of the ratio of the density of the relativistic electrons and of the hot ICM plasma. Right panel: 
Radial profile of the ratio of the pressure of the relativistic electrons and the hot ICM plasma. 



The deprojected densities of the relativistic electrons are 
shown on the right panel of Fig [5] On the left panel of Fig. [6] 
we show the radial distribution of the ratio of the density of 
relativistic electrons and of the hot ICM plasma. On the right 
panel, we show the radial distribution of the ratio of the non- 
thermal to thermal electron pressure, where the pressure of the 
non-thermal electrons was calculated as P = E e /3V. 

The relative density and relative pressure of the relativistic 
electrons peaks strongly in the core of the cluster in the area 
with a radius of 2'. 

3.4. 2D maps of the soft excess emission 

In order to determine the 2D distribution of the soft excess 
emission and to look for its spatial variations, we extract 25 
spectra using XMM-Newton EPIC from a rectangular grid with 
a bin-size of 3' in the 15'xl5' region centred on the core of the 
cluster. For each bin, we compute a spectral redistribution file 
and an ancillary response file. We fit the spectrum of each bin 
individually. We model the hot ICM plasma emission with the 
multi-temperature wdem model and the soft excess emission 
with a thermal component (kT = 0.19 keV, metallicity Z = 0.3 



solar). As shown on the left panel of Fig. [7J we confirm the 
peaked distribution of the emission of the soft component. We 
also find no significant azimuthal variations in the soft excess 
emission. On the central panel of Fig. [7j we show the ratio of 
the emission measure of the warm gas and of the hot ICM. We 
again see that although the soft excess emission strongly peaks 
in the cluster core, its relative contribution is there the smallest 
and increases with the cluster radius. Both the soft excess and 
the cluster emission peak on the central cluster galaxy. On the 
right panel of Fig. [7] we show an optical Digitalised Sky Survey 
image of the region of the cluster used to derive the map of the 
soft excess emission with overplotted extraction regions and 
X-ray contours from Suzaku XIS 1 . 

3.5. Search for line emission from the warm gas with 
RGS 

Both the radial profiles and the 2D maps of the soft excess 
emission show that the emission is strongly peaked on the clus- 
ter core. If the soft excess emission is due to cooling gas in the 
core of the cluster, we could detect its presence via redshifted 
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Fig. 7. Left panel: Map of the emission measure of the soft excess in the inner 15'xl5' region of the cluster. The units are 10 63 
cirr 3 arcmkr 2 . Central panel: Map of the ratio of the hot ICM plasma emission measure and soft excess emission measure in 
the same region as the previous map. Right panel: Optical Digitalised Sky Survey image of the central region of the cluster that 
we used to derive the map of the soft excess emission. The extraction regions for the maps and the X-ray contours from Suzaku 
XIS 1 are overplotted. Both the soft excess and the cluster emission peak on the central cluster galaxy. 
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Fig. 8. RGS spectrum extracted from a 4' wide stripe centred on 
the cluster core. The thick and the thin lines depict the best fit 
models of the RGS data without and with the warm gas compo- 
nent, respectively. Note the mismatch between the two models 
at ~ 23 A, the energy of the redshifted O vn emission. 



O vii lines with the RGS detector on XMM-Newton. We fit the 
RGS spectrum extracted from a 4' wide stripe centred on the 
cluster core with the wdem mode l (for the details of the RGS 
spectroscopy of Sersic 159-03 see Ide Plaa et alJl2006l) . We ob- 
tai n a good fit with pa rameters consistent with those reported 
by de Plaa et al. ( 2006h . The RGS spectrum is shown in Fig. [8] 
with the best fit wdem model depicted by the thick line. By 
fitting a Gaussian at the expected wavelength of the O vn inter- 
combination line at the cluster redshift, we find a 2cr upper limit 
for the Ovu line flux of 1.4 x 10~ 6 counts s _1 cirr 2 arcmkr 2 
(assuming O vn emission from a circular region with a radius 
of 2'). 

Using the sum of the deprojected emission measures of the 
warm gas obtained for the 3 extraction regions within the in- 



ner 2' of the cluster (T e xc2' = 3.7 x 10 65 cm" 3 ), we model 
the thermal spectrum of the soft excess emission. In the model 
we assume the parameters (kT = 0.19, Z = 0.3) used in 
Sect. 13.31 to determine the deprojected radial distribution pro- 
files. The thin line in Fig. [8] shows the best fit model to the 
RGS spectrum with the warm gas component included in the 
model. The model for the warm gas predicts an O vn line flux 
of 3.2 x 10~ 6 counts s cirr 2 arcmin~ 2 . This line flux, and thus 
the presence of warm gas with kr = 0.2 keV, metallicity of 0.3 
solar, and mass of M warm = 4.4 x 10 11 M , can be ruled out 
using the RGS at the 4.6<x confidence level. 



4. Discussion 

4.1. Presence of the soft excess 

Although the systematic uncertainties in the determined soft 
excess flux are large (about 30%), the Suzaku observation con- 
firms the presence of the soft excess emission in Sersic 159- 
03 and the flux derived using Suzaku data is within the sys- 
tematic uncertainties consistent with that found with XMM- 
Newton. The presence of the soft excess in this cluster is now 
established by Suzaku, by two independent observations with 
XMM-Newton and by a previous detection with ROSAT. The 
flux of the soft excess emission in the 9'— 12' radius is at least 
as large as the total flux of the Galactic foreground emission 
(it is more than 5 times larger than the flux of the hard dis- 
tant foreground component) and the soft excess strongly peaks 
on the cluster core. The observed high flux and the radial dis- 
tribution of the excess emission rule out the possibility that 
it is due to a wrong subtraction of Galactic foreground emis- 
sion, o r due to charge exchange emiss ion as previously sug- 
gested ( Bregman & Llovd-Daviesll2006l) . We conclude that the 
soft excess emission is present and is intrinsic to the cluster of 
galaxies Sersic 159-03. 
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4.2. Thermal or non-thermal origin? 
4.2.1. Thermal models 

We find that the soft excess emission peaks strongly on the cen- 
tral cD galaxy of Sersic 159-03. If the observed soft excess 
would be due to warm-hot intergalactic medium (WHIM) in- 
tercluster filaments, then the densest central part of the filament 
would have to be exa ctly along our line of sig ht to the cluster 
core and as shown by lBonamente et al.1 (120051) . for an assumed 
density of n = 10~ 4 cm -3 in its centre, the filament would be 
several hundred megaparsecs long. For smaller assumed den- 
sities it would be even longer than 1 Gpc. This scenario is not 
consistent with current cosmological models and observations. 
By comparing the observed properties of the soft excess emis- 
sion reported in the liter ature with the predicti o ns based on 
model s of the WHIM bv lCen & Ostrikerl d 19991) . iMittaz et al. 
(2004) concluded that the emission of the WHIM in current 
models is 3—4 orders of magnitude too faint to explain the soft 
excess. 



The preferred interpretation of lBonamente et alj d2005l) for 
the soft excess is warm gas with a density of 10~ 3 -10~ 2 cm" 3 
in a clumpy intracluster medium with a volume filling fac - 
tor much smaller than 1 as suggested by ICheng et al.l (2005). 
According to this interpretation the warm gas would have a 
cooling time shorter than the Hubble time and the soft excess 
would be a dynamical phenomenon observable due to recent 
infall of smaller groups. However, only a model where the sub- 
group is falling in just along our line of sight would be consis- 
tent with the centrally peaked distribution of the soft emission, 
without significant azimuthal variations. But in that case, we 
would not observe the soft excess out to large radii as we do. 
This model is not consistent with the observations. 

From observations, the thermal interpretation can only be 
confirmed by the presence of O vn line emission at the cluster 
redshift. However, the Suzaku observation does not confirm the 
presence of these lines. The XMM-Newton RGS data rule out 
the presence of the large amount of warm gas in the cluster core 
deduced from the deprojected radial profiles assuming Z = 0.3 
and kT = 0.19 at the 4.6<x confidence level. Using the RGS, we 
found an upper limit for the O vn line flux in the cluster core of 
3.2 x 10~ 6 counts s~ 



4.2.2. Non-thermal model 

The data are consistent with the non-thermal description of the 
soft excess. If the soft excess is of non-thermal origin, the total 
energy in the relatitivistic electrons producing X-ray emission 
between 0.3 and 10.0 keV within 8' (600 kpc), calculated us- 
ing equation ([TJi, and using the luminosities and photon indeces 
determined in Sect. 13. II is ~2 x 10 60 ergs. If we assume that the 
relativistic electrons have the same differential electron distri- 
bution down to energies with Lorentz factors of y ~ 200, be- 
low which the electrons loose energy rapidly due to Coulomb 
losses, the total energy in relativistic electrons will still not ex- 
ceed 1 x 10 61 ergs. The total thermal energy within the radius 
of 600 kpc of the cluster is 3 x 10 63 ergs, which means that 
even if the energy in the relativistic ions is as much as ~30 



times higher than that in relativistic electrons, the total energy 
in cosmic ray particles will only account for 10% of the thermal 
energy of the ICM. The distribution of the relativistic electrons 
peaks on the central cluster galaxy. Outside the region with a 
radius of 2' (150 kpc), the density of non-thermal electrons rel- 
ative to the electron density of the hot ICM stays constant. The 
relative density and pressure of the relativistic electrons at the 
central cluster galaxy are « 10 times higher than in the outer 



parts. This result is in agreement with that of Bon amente et al 
(20051). 



The high non-thermal electron density in the cluster core 
might have been provided by an AGN in the central cD 
galaxy. Our best fit power-law photon indices in the inner 
and outer extraction region of T = 2.11 and F = 2.41, re- 
spectively, correspond to differential electron number distri- 
butions of ju = 3.22 and fi = 3.82, respectively. These are 
steeper than the distribution of the Galactic cosmic-ray elec- 
trons (// ~ 2.7). The steepening of the power-law distribu- 
tion might indic ate that the r elativistic electrons suffered ra- 
diative losses dSarazin|ll999l) . Relativistic electrons emitting 
in the observed energy range via IC emission have relatively 
long lifetimes. The respective inverse compton tic and syn- 
chrotron lifetimes are: tic = 2.3 x 10 9 (y/10 3 )~ | (1 + z)~ 4 yr a nd 
f syn = 2.4 x 10 10 (y/lO 3 )- 1 (B/lfiG)- 2 yr (e.g. ISarazinlll999l) . 

The relativistic particles could be seeded by jets of active 
galaxies, then transported outwards whilst undergoing in-situ 
second order Fermi a cceleration by turbulent Alfven waves 
dLieu & Ouenbv 2006). The presence of relativistic particles 
at large distances from the cluster core might also be explained 
by shock accelerati on. For examp le by accretion shocks or as 
suggested by Ide Plaa et al. ( 2006), by shocks associated with 
ram pressure stripping of infalling galaxies. Diffuse shocks ac- 
celerate both electrons and ions and the energy in i ons could b e 
~20-50 times larger than the energy in electrons QBellll 19780 . 
If we take this into account, then the pressure in the cluster 
center will be close to equipartition with the hot ICM pressure. 
However, if a significant fraction of the relativistic electrons 
associated with the central peak is originating in AGN elec- 
tron/positron jets then the pressure in the relativistic particles 
might still be below pressure equipartition in the cluster core. 

Recently, iBagchi et al. I d2006h reported that the ring-shaped 
non-thermal radio emitting structure found at the outskirts of 
the rich cluster of galaxies Abell 3376 may trace the shock 
waves of cosmological large scale matter flows. The giant radio 
structures in Abell 3376 are probably due to merger or accre- 
tion shocks near the virial infall region of the cluster, which 
is the transition zone between the hot cluster medium and the 
WHIM. Radio sources, like that in Abell 3376 may also be 
the s ites where magneti c shocks accelerate cosmic ray parti- 
cles (Bagchi et al. 2006}). The example of Abell 3376 and of 
other clusters with radio halos shows that in clusters of galaxies 
electrons are accelerated to relativistic energies with a Lorentz 
factors of the order of 10 4 even at large radii. The strong soft 
excess in the outer parts of Sersic 159-03 might be due to accre- 
tion shocks. In order to get a more complete picture about the 
physical processes in Sersic 159-03, it is important to obtain 
radio data of the cluster, that will probe the population of the 
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relativistic electrons with higher energies than those the pres- 
ence of which we infer from the detected soft excess emission. 

Clusters of galaxies are more complicated than previously 
thought. If the soft excess emission is of non-thermal origin, 
as our data suggest, then the non-thermal component needs to 
be taken into account in the overall mass and energy budget 
of galaxy clusters. Reduction of the emission measure of the 
hot cluster gas in the outer parts of the cluster by 40% means 
that the amount of gas in the outskirts is 20% less than indi- 
cated by models which do not take into account the presence of 
the non-thermal emission. The cluster masses are calculated by 
assuming hydrostatic equilibrium, not taking into account the 
possible non-thermal pressure. The non-thermal pressure may 
be high enough to affect the estimates of the total cluster mass. 

If a substantial fraction of the central 2-10 keV emission 
of some clusters of galaxies is of non-thermal origin, it will 
not only account for the detection of soft and hard excesses, 
but also for the significantly less than expected Sunyaev- 
Zel'dovich effect, that emerged from a comparison of WMAP 



and X-ray data for a large s ample of nearby clusters (ILieu et al 



2006; Afshordi et al. 2006). Because for an equal pressure 
of thermal and relativistic electrons the Sunyaev-Zel'dovich 
decrement is much sm aller due to the relativistic electrons (e.g. 
EnBlin & Kaisei1l2000l) . 



5. Conclusions 

We have analysed new Suzaku XIS1 spectra and archival 
XMM-Newton EPIC and RGS spectra of the cluster of galaxies 
Sersic 1 59-03 . We found that: 

- The Suzaku observation confirms the presence of the soft 
excess emission and its derived flux is consistent with the 
values determined using XMM-Newton. Suzaku does not 
confirm the presence of the redshifted O vn lines in the clus- 
ter. 

- Radial profiles and 2D maps show that the soft excess emis- 
sion has a strong peak at the position of the central cD 
galaxy and the maps do not show any significant azimuthal 
variations. Although the soft excess emission can be fitted 
equally well with both thermal and non-thermal models, 
the spatial distribution of the soft emission is not consistent 
with any of the viable thermal models: intercluster WHIM 
filaments, or models of clumpy warm intracluster gas asso- 
ciated with infalling groups. 

- Using the XMM-Newton RGS data we do not confirm the 
presence of the warm gas in the cluster centre as inferred 
from the deprojected EPIC density profiles and we put up- 
per limits on the O vn line emission from the cluster core. 

- The observed properties of the soft excess emission are 
consistent with the non-thermal interpretation. While the 
high density of relativistic electrons associated with the 
peak of the soft emission in the cluster centre might have 
been provided by an AGN in the central cD galaxy, the un- 
derlying population might have been accelerated in diffuse 
shocks. 
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Appendix A: The contamination layer for Suzaku 
XIS1 at the time of the observation 

We analysed the data obtained during a simultaneous XMM- 
Newton EPIC/pn, Chandra LETGS/HRC and Suzaku XIS 1 ob- 
servation of the blazar PKS 2155-304. The observation was 
performed on May 1-2, 2006, only 3 days after our observa- 
tion of the cluster of galaxies Sersic 159-03. This observation 
allows us to determine the thickness of the contaminating layer 
in the centre of the field of view of Suzaku XIS 1 at the time 
of our observation. Combining the best fit central value with 
the radial contamination profile, that is determined by the XIS 
team, we can determine the contaminating column in each anal- 
ysed region at the time of our observation. 

A.1. Data reduction 

We used version 1.2 cleaned event files to extract the XIS1 
spectrum. The spectrum was extracted from a circular region 
with a radius of 250 detector pixels. Although the region is not 
fully covered in the observation due to the 1/4 window obser- 
vation mode, the loss of flux is only a few percent and inde- 
pendent of the energy. Hence, we can neglect this effect. Since 
PKS 2155-304 is a highly variable object, we only used the 
time intervals when both XMM-Newton EPIC/pn and XIS1, or 
Chandra LETGS and XIS1 data were available. The RMF file 
for XIS 1 was created using "xisrmfgen". We used the ARF file 
"ae_xil_xisnom6_20060615.arf" (this file is distributed by the 
XIS team), which does not take into account the contamination 
on the optical blocking filter (OBF). 

The XMM-Newton data were reduced using SAS 7.0.0. 
After selecting only those time intervals when both Suzaku and 
XMM-Newton were observing the source, we analysed 8.6 ks 
of data from both satellites. The EPIC/pn observation was per- 
formed in the small window mode. We verified that the ob- 
tained spectrum is not affected by pile-up. 



Chandra LETGS/HRC data were reprocessed according to 
the standard procedure of CIAO 3.3.0 with CALDB 3.2.3. 
After selecting only those periods when both LETGS and XIS 1 
data were available, we analysed 17.8 ks of data from both 
satellites. The spectra of plus and minus orders were merged. 
We created ±lst— 8th order response files (RMFs and Grating 
ARFs). In order to take into account the overlap of the dis- 
persed photons of each order, we adopted the sum of ±lst-8th 
order response files as the response for the LETGS spectrum. 

A.2. Simultaneous Suzaku XIS1 and XMM-Newton 
EPIC/pn observation 

We fitted the spectra using the SPEX spectral fitting package. 
The source spectrum was fitted with a broken power-law. For 
the Galactic absorption model, we used the "hot" model of 
SPEX, where we set the abundances in the absorber to the 
proto-solar values of Lodders (2003). We fixed the hydrogen 
column density in the model to the best determined value of 
Nu = (1.36 + 0.10) x 10 20 cm 4 , which wa s found using a ded- 



icated H i observation toward this blazar dLockman & Savage 
1995b . We fitted the EPIC/pn spectrum in the 0.35-10 keV 
band. Our best fit low energy power-law photon index is T — 
2.61 + 0.01. The break in the spectrum is at 2.0 + 0.3 keV and 
its best fit value is AF = -0.1 1 + 0.03. 

We used the best fit model determined using EPIC/pn for 
the XIS1 data. The XIS1 spectra were fitted in the 0.35-7 keV 
band. In order to find the best fit absorption column on the 
OBF of XIS 1 , we included in the model free C i and O i ab- 
sorbing columns. Their best-fit values are Nc = (4. 1 1 ± 0.06) x 



-2 



10 1S cm- z and N = (6.2 + 0.5) x 10 1 ' cm 

The contamination values determined using the EPIC/pn 
are only slightly lower than the nominal values from the stan- 
dard contamination file for September 1 . 



A.3. Simultaneous XIS1 and Chandra LETGS/HRC 
observation 

We fitted the LETGS and the XIS 1 spectra simultaneously us- 
ing the SHERPA spectral fitting package. The fitting range was 
0.17-2.0 keV for LETGS and 0.35-7.0 keV for XIS1. The 
spectrum of PKS 2155-304 was modelled as a power-law con- 
volved with Galactic absorption (wabs model). The model was 
additionally convolved with C i and O i absorption by the con- 
tamination on the OBF of XIS 1 . Since the most probable com- 
position of the contaminant is C24H38O4, the C/O ratio was 
fixed in the fitting process to 6. We fixed the hydrogen col- 
umn density in the model to the best determined Galactic value 
of A^h = (1.36 + 0.10) x 10 2() cm- 2 . The spectra were well 
fitted with a photon index of F = 2.54 + 0.01. The best fit con- 
taminating columns are Nc = (3.78 + 0.05) x 10 18 cirT 2 and 
No = 6.30 x 10 17 cm" 2 Note that the photon index and the 
flux are consistent between LETGS and XIS 1 if they are fit- 
ted separately. We verified that we obtain consistent values for 
the absorbing column also with the A^h as a free parameter. We 
also tried broken power-law model which did not improve the 
fit significantly. 
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The contaminating columns determined using the LETGS 
are very similar to the nominal values from the standard con- 
tamination file for July 1 . 



